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ABSTRACT 



Aims. We study the impact of internal gravity waves (IGW), meridional circulation, shear turbulence, and stellar contraction on the 

internal rotation profile and surface velocity evolution of solar metallicity low-mass pre-main sequence stars. 

Methods. We compute a grid of rotating stellar evolution models with masses between 0.6 and 2.0 Mg taking these processes into 

account for the transport of angular momentum, as soon as the radiative core appears and assuming no more disk-locking from that 

moment on. IGW generation along the PMS is computed taking Reynolds-stress and buoyancy into account in the bulk of the stellar 

convective envelope and convective core (when present). Redistribution of angular momentum within the radiative layers accounts 

for damping of prograde and retrograde IGW by thermal diffusivity and viscosity in corotation resonance. 

Results. Over the whole mass range considered, IGW are found to be efficiently generated by the convective envelope and to slow 

down the stellar core early on the PMS. In stars more massive than ~ 1.6 Mq, IGW produced by the convective core also contribute 

to angular momentum redistribution close to the ZAMS. 

Conclusions. Overall, IGW are found to significantly change the internal rotation profile of PMS low-mass stars. 
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1. Introduction 

The evolution of the surface rotation of low-mass stars along 
the pre-main sequence (hereafter PMS) follows a specific path 
as shown by the data, both rotation periods and v sin i measure- 
ments, collected in young stellar clusters (e.g. Irwin & Bouvier 
2009, and references therein for a review). The rotational prop- 
erties of young stars appear to result from an intricate interplay 
between several physical processes that affect the angular mo- 
mentum gains, losses, and redistribution as the stars evolve along 
the PMS towards the zero age main sequence (hereafter ZAMS). 
These mechanisms can be roughly divided into two classes. The 
first ones result from the connection of the stars to their environ- 
ment (magnetic and dynamic coupling to a circumstellar disk, 
accretion, jets, stellar and disk winds, etc.), and are particularly 
crucial during the T-Tauri phase, when star-disk interaction is 
observed and expected to be strong (see, for instance, Shu et al. 
1994; Matt & Ptidntz 2005; Zanni & Ferreira 2012). The broad 
variety of possible star-environment configurations may, in par- 
ticular, explain part of the large dispersion in rotation rates of 
solar-type stars observed along the PMS and at the arrival on the 
ZAMS (e.g., Stauffer et al. 1985; Irwin & Bouvier 2009). 

The second class is related to stellar secular evolution and 
consists of the (magneto-) hydrodynamical transport mecha- 
nisms that contribute to redistributing angular momentum in- 
side the stars themselves. In the present paper we focus on 
exploring these internal mechanisms once the disk has dissi- 
pated and the accretion process is over, which occurs after 3-10 
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Myr (e.g., Haisch et al. 2001; Hartmann 2005; Hernandez et al. 
2008), i.e., roughly at the time when a radiative core appears 
in the contracting PMS stars. Our aim is to evaluate, in particu- 
lar and for the first time, the interplay between internal gravity 
waves (hereafter IGW), meridional circulation, turbulent shear, 
and stellar contraction during the PMS, considering that IGW 
are one of the best candidate mechanisms to explain the flat an- 
gular velocity profile inside the Sun as revealed by helioseismol- 
ogy (Charbonnel & Talon 2005). This work is also motivated by 
the results of Talon & Charbonnel (2008, hereafter TC08), who 
showed that IGW are efficiently excited inside intermediate- 
mass PMS stars and suggested that waves should efficiently 
transport angular momentum during the PMS evolution, which 
should affect the angular velocity profile at this phase and at the 
arrival on the ZAMS. 



In §2 we introduce the formalism for IGW excitation and for 
the transport of angular momentum through the various mecha- 
nisms considered. We describe in §3 the basic assumptions for 
the present grid of PMS models for low-mass (0.6 to 2 Mq), 
solar-metallicity stars. In §4 we examine IGW generation by 
the convective envelope and the convective core (when present) 
along the PMS evolution for the whole mass range covered by 
the grid. In § 5 we describe the impact of the various interacting 
transport mechanisms on the evolution of the internal rotation 
profile for the various stellar masses, and in § 6 we briefly dis- 
cuss their influence on the surface rotation velocity and lithium 
abundance, as well as on global stellar properties. Conclusions 
are presented in § 7. 
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2. Formalism 

We follow Talon & Charbonnel (2005, hereafter TC05) for the 
treatment of both the excitation of IGW and the transport of 
angular momentum and chemicals by waves, meridional circu- 
lation, and shear turbulence in hydrodynamical stellar models. 
We, however, underline three main improvements over TC05 
paper. First, we consider IGW generated both by the external 
and central convective regions (when present), while only those 
excited by the convective envelope were considered in our pre- 
vious work. Second, the important variations in the stellar struc- 
ture and of the IGW properties along the pre-main sequence (see 
§ 4 and TC08) require that we then compute the wave spectra at 
each evolution time step, while the main sequence computation 
presented in TC05 was based on the wave spectrum of the stellar 
model on the ZAMS. Finally, we account here for both prograde 
and retrograde waves in the whole radiative interior, while only 
the latter ones were considered previously. We recall below the 
formalism (i.e., relevant equations and assumptions) that is in- 
cluded in the evolution code STAREVOL (see TC05, TC08, and 
Mathisetal. 2013). 



2.1. IGW generation 

In this exploratory work we apply the Goldreich et al. (1994) 
formalism as adapted by Kumar & Quataert (1997) to calculate 
the spectrum of IGW excited by Reynolds stress and buoyancy 
in the bulk of convective regions (see e.g. Zahn el al. 1997). We 
do not consider the possible effects of IGWs generated by con- 
vective overshooting plumes, since no analytical prescription is 
available to describe this excitation mechanism (see details and 
discussion in TC05). However as we see in § 2.3.1, we multi- 
ply the IGW flux by a factor 2 in the transport equation in order 
to account for the recent results by Lecoanet & Quataert (2012). 
We treat the waves by assuming that they are pure gravity waves 
(i.e., not modified by the Coriolis acceleration) that only feel the 
entrainment by diff'erential rotation. 

The kinetic energy flux per unit frequency is 



rE{(,co) 
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(1) 



with ^r and [/'(/' + 1)]'^^^/, the radial and horizontal displace- 
ment wave functions normalized to unit energy flux at the edge 
of the considered convection zone, Vc the convective velocity, 
L - aMLjHp the radial size of an energy bearing turbulent 
eddy, Tl ^ L/v^ the characteristic convective time, Hp the pres- 
sure scale height P/pg, and /2„ the radial size of the largest 
eddy at depth r with characteristic frequency of co or higher 
(h^j = Lmin{l,(2a;Ti)"''^^)). The radial and horizontal wave 
numbers (respectively kr and k^) are related by 



e-l"4'^Vl-K 



where A^^ is the Brunt- Vaisala frequency' 
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' The Brunt- Vaisala - or buoyancy - frequency N is given by A'" = 
^r + ^i = tI: (<5(Vad - V) + <^V^), with 6 = {d\npld\nT)p,, ^ = 
(51np/51nyu)p7-, V the logarithmic temperature gradient, and V^ the 
mean molecular weight gradient. 



At the considered convective edge (located at radius rf.,), the 
mean flux of angular momentum carried by a monochromatic 
wave of spherical order £ and local (i.e., emission) frequency co 
{m being the azimutal order), i.e., the momentum flux per unit 
frequency, is related to the kinetic energy flux by 



2m 

0) 



(3) 



(Goldreich & Nicholson 1989; Zahn et al. 1997). The so-called 
angular momentum luminosity at the considered convective edge 
(envelope or core) is obtained after horizontal integration 



■Cj(„,{r^z) ^Anr^JFj.- 



(4) 



2.2. IGW damping 



Deposition of angular momentum (positive or negative) 
within the radiative layers occurs at the depth where in- 
dividual monochromatic waves are eventually damped by 
thermal diffusivity and viscosity in corotation resonance 
(Goldreich & Nicholson 1989; Schatzman 1993; Zahn et al. 
1997). In the present study the local momentum luminosity at a 
given radius r within the radiative region accounts for prograde 
and retrograde waves (i.e., with respectively positive and nega- 
tive m values) generated by both the convective envelope and the 
convective core (if present), i.e.. 



-CAr) = £y,env + -C. 
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where each component is given by 

^ -Cjf,m (i-cz) exp [-T(r, 0-, £)] . 



(5) 
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where 'cz' refers to the interface between the radiative region 
and the corresponding convection zone (i.e., envelope or core). 
The local damping rate 
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(7) 
takes the mean molecular weight stratification into account 
(Goldreich & Nicholson 1989; Schatzman 1993; Zahnetal. 
1997), as well as the thermal and the (vertical) turbulent viscos- 
ity (Kt and v,. respectively). Here, cr is the local Doppler-shifted 
frequency 

o-(r) = w - m [Q(r) - Qez] (8) 

with CO the wave frequency in the reference frame of the corre- 
sponding emitting convection zone that rotates with the angular 
velocity ficz- 

As can be seen from these expressions, angular momen- 
tum redistribution by IGW within the radiative region is dom- 
inated by low-frequency (cr «: N), low-degree waves; indeed, 
those penetrate deeper, and their prograde and retrograde com- 
ponents experience strong differential damping, as required to 
produce a net momentum deposition. In contrast, high-degree 
waves are damped closer to the convection zone (since damping 
oc [£(£ + 1)]2), and high-frequency waves experience less differ- 
ential damping. 
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2.3. Global transport of angular momentum by IGW, 
meridional circulation, and shear turbulence 

2.3.1. General equations 

We assume solid-body rotation in the convective regions. In the 
stellar radiative regions, the evolution of angular momentum 
through advection by meridional circulation, diffusion by shear 
turbulence, and deposit or extraction by IGW follows the general 
expression below (e.g., Talon & Zahn 1998): 
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where U is the radial meridional circulation velocity, v,, the tur- 
bulent viscosity due to differential rotation, and p the density. 
We have added a factor 2 in the last term to account for the 
study by Lecoanet & Quataert (20 1 2), who predict the IGW flux 
due to turbulent convection to be a few to five times larger than 
in previous estimates by, e.g., Goldreich & Kumar (1990) and 
Goldreich et al. (1994). However as we see in § 5.1.2, our con- 
clusions are not sensitive to this multiplication factor. 

Following Decressin et al. (2009) and Mathis et al. (2013) 
we integrate Eq. 9 over an isobar enclosing the mass m (r) to 
obtain the expression of the total flux (loss or gain) of angular 
momentum carried by the considered transport processes: 
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where the fluxes driven by meridional circulation, vertical shear- 
induced turbulence, and IGWs are, respectively. 
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For the treatment of horizontal turbulent viscosity, we follow 
Zahn (1992), again as in TC05: 



Dh = Vft = — 



1 d(pr^U) Udlnr^D. 



3pr dr 
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withC/, - 1. 

The influence of the prescriptions assumed for Di. and D/, 
will be investigated in a future paper. 

2.4. Transport of chemicals 

We treat the transport of chemical species in the radiative regions 
as a diff'usive process through the combined action of meridional 
circulation and shear-induced turbulence (Chaboyer & Zahn 
1992). The effective diffusion coefficient is written 



£>eff = 



30 Dh 



(15) 



where £);, is the horizontal component of the turbulent diffiisivity 
(see Eq. 14). 

In the present study we neglect atomic diffusion, whose ef- 
fects require much longer timescales to develop compared to the 
very short duration of the pre-main sequence phase. We also ne- 
glect possible wave-induced turbulence. Therefore the expres- 
sion for the transport of chemicals (here, the mass fraction X of 
the element in the stellar radiative region writes as (see e.g. 
Meynet & Maeder 2000): 
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oM,- 
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where dMr - A-npr^dr, and the last term accounts for nuclear 
destruction or production of the considered element. 
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2.3.2. Meridional circulation 



As can be seen in Eq. 9 the transport of angular momentum 
through meridional circulation is treated as an advective process. 
As in our previous studies we apply the formaUsm developed by 
Zahn (1992), Maeder & Zahn (1998) and Mathis & Zahn (2004, 
see also Decressin et al. 2009). 

2.3.3. Shear-induced turbulence 

Shear-induced turbulence is assumed to be highly anisotropic. 
Following TC05 we assume that the turbulent diffusion coeffi- 
cient equals turbulent viscosity and use the corresponding ex- 
pression given by Talon & Zahn (1997), i.e.. 



A 



8 Ri„it(rdn/dry 

5Nll{K + Dh) + NllD,, 



(13) 



that considers the weakening effect of thermal diffusivity (Kj) 
on the thermal stratification and of horizontal turbulence (D/,, 
see below) on both the thermal and mean molecular weight strat- 
ifications. 



3. Stellar models 

3.1. Input physics and basic assumptions 

We focus on the pre-main sequence evolution of solar- 
metallicity stars in the mass range between 0.6 and 2.0 Mq. We 
adopt the solar composition of Asplund et al. (2009). Opacity 
tables are updated accordingly both at high and low tem- 
perature respectively from OPAL and Wichita websites^ (see 
e.g. Iglesias & Rogers 1996; Ferguson et al. 2005). The mixing 
length parameter qtmlt = 1 -63 is calibrated so that our standard 
(i.e., non rotating) 1 Mq, Zq model fits the solar radius, effective 
temperature, and luminosity at the age of the sun. Convection 
zone bounderies are defined by the Schwarzschild criterion, and 
we do not account for convective overshoot. 

Computations are performed with the stellar evolution code 
STAREVOL (see e.g. TC05, Lagarde et al. 2012). Initial models 
are totally convective polytropic stars, with central temperature 
lower than 10'' K (i.e., deuterium burning has not yet occurred). 
We follow the PMS evolution along the Hayashi track up to the 
arrival on the ZAMS that we define as the point where the ratio 
between central and surface hydrogen abundance reaches 0.998. 
The stellar mass is assumed to be constant during that phase (i.e., 
no accretion nor mass loss). For each stellar mass we compute 



^ http : //adg . llnl . gov/Research/OPAL/opal . html ; 
http : //webs .Wichita . edu/physics/opacity/ 
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Fig. 1. PMS tracks in the Hertzsprung-Russel diagram for solar metallicity stars with initial masses between 0.6 and 2.0 Mq (classical 
models are shown here) and properties of the external convective layers. Colors indicate the radial extent of the convective envelope 
(top left panel), the temperature at its bottom (top right panel), the maximal convective flux (bottom left panel), and the thermal 
diffusivity below the envelope (bottom right panel). Dashed lines connect points with similar values for these quantities, and the 
colored axes are in cgs units. The dotted parts of the tracks correspond to the phase when the stars are still fully convective 



classical models (i.e., without any transport of angular momen- 
tum nor of chemicals) as well as rotating models with and with- 
out IGW. We neglect the hydrostatic effects of the centrifugal 
force in all our rotating models but two; we discuss the impact 
of this simplification in § 5. The evolution tracks of the classical 
models in the Hertzsprung-Russel diagram are shown in Fig. 1. 



3.2. Initial internal and surface rotation 

We assume solid-body rotation while stars are fully convective 
(which corresponds to the dotted part of the tracks in Fig. 1) 
and we start computing the evolution of surface and internal 
rotation under the action of stellar contraction, meridional cir- 
culation, turbulence, and IGW when the radiative core appears, 
which happens at ages between ~ 0.5 and 7.5 Myr for the mass 
range considered (see T(core) in Table 1), and at ~2.5 Myr for 
the 1.0 Mq model. At that time, most or even all low-mass 
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Fig. 2. Same as Fig. 1, but for the properties of the convective core. Here dotted parts on the tracks indicate the phase when the 
convective core is not yet present 



stars have already lost their disks as shown by observations in 
very young clusters (e.g. Haisch et al. 2001; Hartmann 2005; 
Hernandez et al. 2008). For all stellar masses, we choose the ini- 
tial rotation velocity at the moment when the radiative zone ap- 
pears to be equal to 5% of the critical velocity of the correspond- 
ing model (Vert — a/I^' ^^^ Table 1). This corresponds ap- 
proximately to the median of the observed distribution in young 
open clusters (see Fig. 13 and § 6 for discussion). We assume 
that there is no more coupling between the star and a potential 
disk beyond that evolution point. The surface of the star is then 
free to spin up, and we do not apply any magnetic braking. The 



influence of the initial rotation velocity, of the disk lifetime that 
affects the moment when a PMS star starts spinning up, as well 
as that of magnetic wind braking that may affect the rotation rate 
at the arrival on the ZAMS, will be investigated in a forthcoming 
paper 

4. IGW generation along the PMS evolution for all 
grid models 

The internal structure strongly changes as low-mass stars evolve 
along the PMS. This implies strong variations in the quantities 
that are relevant to IGW generation and momentum transport, as 
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Table 1. Properties of the different models computed without rotation (std), with rotation but without IGW (rot), and with rotation 
taking into account IGW (igw); for the 1.0 Mq star a couple of models were computed taking into account the hydrostatic effects of 
rotation (rot+hydro and igw+hydro). For each model we give: lifetime on the PMS, age at which the radiative core appears, initial 
rotation velocity, surface rotation rate and rotation period when the radiative core appears (taken at 5% of critical rotation velocity of 
the corresponding model), and surface rotation velocity, surface rotation rate, rotation period, surface lithium abundance, effective 
temperature, and luminosity at the arrival on the ZAMS 



Star 




T(pms) 


T(core) 


V.,„r/ 


fi™,-/ / Ho 


l^init 


v.„„-/- 


n™,f/no 


PzAMS 


N(Li) 


Teff' 


log L/Lo 


Mo 




Myr 


Myr 


init, km.s"' 


mit 


days 


zams, km.s"' 


zams 


days 


zams 


zams, K 


zams 


0.6 


std 


191 


7.47 


- 


- 


- 


- 


- 


- 


-4.96 


4118 


-1.14 




rot 


194 


7.47 


14.5 


6.60 


3.73 


36.3 


34.5 


7.30 


-5.08 


4127 


-1.14 




igw 


217 


7.47 


14.5 


6.60 


3.73 


37.3 


35.6 


7.12 


-5.40 


4121 


-1.14 


0.9 


std 


75 


3.00 


- 


- 


- 


- 


- 


- 


2.66 


5246 


-0.366 




rot 


76 


3.00 


14.1 


4.16 


5.99 


48.5 


30.7 


0.83 


2.64 


5246 


-0.366 




igw 


89 


3.00 


14.1 


4.16 


5.99 


55.7 


35.2 


0.72 


2.64 


5246 


-0.366 


1.0 


std 


57.5 


2.43 


- 


- 


- 


- 


- 


- 


2.96 


5608 


-0.152 




rot 


53.9 


2.43 


14.2 


3.76 


6.84 


53.1 


29.8 


0.76 


3.02 


5608 


-0.152 




rot+hydro 


59.4 


2.43 


14.2 


3.76 


6.84 


52.9 


29.7 


0.76 


3.12 


5589 


-0.152 




igw 


59.5 


2.43 


14.2 


3.76 


6.84 


63.1 


35.5 


0.68 


2.95 


5608 


-0.152 




igw+hydro 


59.2 


2.43 


14.2 


3.76 


6.84 


63.1 


35.3 


0.72 


2.95 


5584 


-0.156 


1.2 


std 


38.6 


1.58 


- 


- 


- 


- 


- 


- 


3.17 


6217 


0.235 




igw 


43.5 


1.58 


17.2 


2.99 


8.03 


86.7 


38.6 


0.86 


3.17 


6217 


0.235 


1.4 


std 


26.5 


1.09 


- 


- 


- 


- 


- 


- 


3.22 


6828 


0.574 




igw 


29.4 


1.09 


13.6 


4.08 


6.02 


109 


39.9 


0.64 


3.22 


6838 


0.572 


1.6 


std 


19.6 


0.81 


- 


- 


- 


- 


- 


- 


3.23 


7757 


0.835 




rot 


19.3 


0.81 


13.0 


1.85 


13.7 


100 


39.2 


0.71 


3.23 


7757 


0.837 




igw 


21.2 


0.81 


13.0 


1.85 


13.7 


120 


42.2 


0.57 


3.23 


7783 


0.829 


1.8 


std 


15.0 


0.62 


- 


- 


- 


- 


- 


- 


3.24 


8666 


1.04 




rot 


15.0 


0.62 


12.4 


1.43 


16.9 


112 


37.6 


0.68 


3.24 


8666 


1.04 




igw 


16.1 


0.62 


12.4 


1.43 


16.9 


128 


43.9 


0.57 


3.24 


8676 


1.04 


2.0 


std 


11.8 


0.48 


- 


- 


- 


- 


- 


- 


3.24 


9479 


1.23 




rot 


11.7 


0.48 


20 


3.32 


7.62 


124 


40.7 


0.62 


3.24 


9473 


1.23 




igw 


12.1 


0.48 


20 


3.32 


7.62 


139 


46.0 


0.54 


3.24 


9492 


1.23 



depicted in Figs. 1 for the properties of the convective envelope 
and 2 for the core. All quantities are given in cgs units. 

Stars are first fully convective and a radiative core appears 
along the Hayashi track as they contract and heat (Fig. 1). The 
thickness of the convective envelope decreases, and the temper- 
ature at its base increases as the stars move towards higher effec- 
tive temperatures. Due to central CNO-burning ignition on the 
final approach towards the ZAMS a convective core develops 
(Fig.2). 

One can follow the evolution along the tracks of the maxi- 
mum convective flux {Fc = Cpp\\AT) inside the external and 
central convective regions, which directly affects the energy flux 
associated to a given frequency (see Eq. 1). Wave excitation is 
stronger when the convective length scale {tc - Inrcz/aHp) 
is larger, but decreases when the turnover timescale (r^ - 
OMuHplvc) becomes too large. The combination of these two 
factors induces large differences in the overall efficiency of wave 
generation as the internal structure evolves. This is well illus- 
trated in Fig. 3 that shows the luminosity spectrum of IGW gen- 
erated by the external convection zone in the 1 M© model at four 
ages on the PMS. One sees clearly that wave-induced transport 
is dominated by low-frequency waves (i.e., < 3.5 yuHz). High 
degree waves at low frequencies do not contribute much to the 
transport of angular momentum even though their excitation flux 
is important : indeed they are essentially damped near the con- 
vective envelope edge. In Fig. 4 colors along the tracks indicate 
the net momentum luminosity Xy (see Eq. 4) of IGWs generated 
by the external and internal convective regions. 

In the case of external convection, the net momentum lu- 
minosity -Cisiirf rapidly increases as the excitation of IGW 
strengthens up when stars evolve towards higher effective tem- 



a39 



,2,2 



perature, and reaches maximum values as high as 10 g.cm^s^ 
around Teff ~ 6200 K. Stars with initial masses lower than 1.3 
Mq never reach this effective temperature and the correspond- 
ing X.j.siirf remains always below this maximum and shows only 
a monotonic increase along the PMS. On the other hand in the 
more massive models the convective envelope keeps shrinking 
in size and Xy..!»r/ decreases when Tp/y increases above 6200 K. 
This behavior confirms TC08 findings for intermediate-mass 
PMS stars, and is very similar to the £.j plateau we found for 
Pop I and Pop II main sequence stars (Talon & Charbonnel 2003, 
2004), which share very similar convective properties with PMS 
stars in the same T^ff range. 



IGW are also emitted from the convective core at the end of 
the PMS. The more massive the star, the more the convective 
core expands, and the stronger the corresponding wave excita- 
tion. We note however from Fig. 4 that wave excitation by the 
convective core (when present) is generally much less efficient 
than that of the convective envelope. The ratio between -Cj,core 
and X.j,surf is shown in Fig. 5 as a function of T(.yy for the various 
models. For stars with masses below 1.4 Mq, Jlj^core is always ~ 
5-6 order of magnitude lower than Hj^smf- These two quanti- 
ties reach similar orders of magnitude only very close from the 
ZAMS for stars more massive than 1.6 Mq. Therefore and as we 
shall see below, the impact of IGW on the internal rotation pro- 
file along the PMS will be dominated by the waves emitted by 
the convective envelope. 
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Fig. 3. Angular momentum luminosity at the base of the convective envelope of IGW generated by Reynolds-stress in the external 
convective layers as a function of emission frequency o) and degree L The color axis is in log and white areas correspond to 
logXy,.s«r/ < 22. The plots are shown for the 1 M© model at four ages along the PMS (5.8, 14, 35, and 55 Myr from top left to 
bottom right; the corresponding values of T^ff are 4277, 4357, 5560, and 5612 K). 



5. Evolution of the internal rotation profile 

5.1. The case of the 1 M^ star 

Figure 6 depicts the evolution along the PMS of the rotation pro- 
file inside the 1 M© star for two cases: when angular momentum 
transport is operated solely by meridional circulation and shear 
turbulence (bottom panels), and when angular momentum depo- 
sition by internal gravity waves is taken into account in addition 
to the hydrodynamic processes (top panels); the rotation profile 



is shown at different ages as a function of both relative mass 
fraction and reduced radius (left and right panels respectively). 
The decomposition of the total flux of angular momentum into 
the various components driven by meridional circulation, shear 
turbulence, and IGW (when accounted for; see Eqs. 10, 11, and 
12 respectively) is shown in Fig. 7 at three ages along the PMS. 
Meridional circulation currents are shown at the same ages in 
Fig. 8; clockwise currents (matter flowing from the equator to 
the pole and resulting in deposition of angular momentum in- 
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Fig. 4. Same as Fig. 1, but with colors indicating the total momentum flux carried by IGW generated by the convective envelope 
(left) and the convective core (right). Dotted lines indicate the region where the stars are fully convective (left) or have no convective 
core (right) so that no IGW can be generated. In the left panel the vertical dashed Unes connect models where the excitation has the 
same value: log mTAoj, I, m)\) = 28, 30, 32, 34, 36, 38. 
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Fig. 5. Ratio of the total momentum luminosity carried by IGW 
generated by the convective core and the convective envelope for 
the PMS models of various masses 



wards) and counterclockwise ones (carrying angular momentum 
outwards) are drawn in blue and red respectively. 



5.1 .1 . Transport of angular momentum by meridional 
circulation and shear turbulence only 

When only meridional circulation and shear turbulence are ac- 
counted for, differential rotation rapidly develops inside the ra- 
diative region as the surface rotation velocity increases due to 
stellar contraction (Fig. 6, bottom plots). This behavior as well 
as the rotation profile at the aiTival on the ZAMS are similar to 
the results of Eggenberger et al. (2012) for their rotating 1 M© 
model computed with similar assumptions. 

As can be seen in Fig. 7 (right plots) for this model with- 
out IGW, the transport of angular momentum is dominated by 
meridional circulation all along the PMS, while the contribution 
of shear turbulence is negligible (the flux of angular momen- 
tum by turbulence Fy is indeed ~ 2 orders of magnitude lower 
than the flux driven by meridional circulation Fmc)- The number 
of circulation loops evolves with time (Fig. 8, bottom plots; see 
also Fig. 7): In the early stages (14.9 Myrs, left panel), the cir- 
culation consists of a single counterclockwise current that trans- 
ports matter inward along the rotational axis and outward in the 
equatorial plane; later on (33.5 Myrs, middle panel) a clockwise 
loop appears in the central regions; finally an additional coun- 
terclockwise loop shows up when the convective core develops 
(55 Myrs, right panel). 



5.1 .2. Impact of internal gravity waves 

The evolution of the internal rotation profile changes drastically 
when IGW are taken into account in conjunction with meridional 
circulation and shear turbulence, as can be seen in Fig. 6 (top 
plots). As akeady discussed in § 4, the emitted wave spectrum 
strongly evolves with the stellar structure along the PMS. IGW 
are first emitted by the receding convective envelope, and much 
later by the convective core when it appears during the final ap- 
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Fig. 6. Evolution along the PMS of the rotation profile in the 1 Mq models computed with and without IGW (top and bottom 
respectively) as a function of relative mass fraction and radius in solar units (left and right respectively). The curves are labeled 
according to age. On each plot the left and right scales give Q in //Hz and in solar units respectively 



proach towards the zams. In the case of the 1 Mq model, IGW 
emitted by the convective core play actually no role since their 
luminosity is extremely low (see Fig. 5 and discussion in § 4). 
Therefore the following discussion refers only to those emitted 
by the envelope. 

In order to understand wave-induced transport, we must also 
focus on the important quantities for wave damping in the radia- 
tive layers, namely the Brunt-Vaisala frequency A^- and the ther- 
mal diffusivity Kj'- For a given differential rotation within the 
radiative layers, low-frequency (i.e., with oj < 3.5 //Hz) and/or 
large degree waves that dominate the angular momentum trans- 
port are damped very efficiently close to the convective edges 
when Nj is too small or when Kj is too large (see Eq. 7). Fig. 9 
and 14 show the radial profiles of these two quantities in the ra- 
diative layers of the 1 Mq model at various ages (see also Fig. 1 
and 2 that show the variations along the evolution track of the 
value of Kt just below the convective envelope and above the 



convective core). At all ages A^^ drops near the stellar center and 
the convective edges; in addition its value at a given depth in- 
creases with time along the PMS as a result of the stellar con- 
traction that leads to an increase of gravity and a decrease of the 
pressure scale height as the star evolves. On the other hand the 
value of Kt just below the convective envelope also increases as 
the star contracts and move towards higher effective temperature; 
this implies stronger damping of all the waves (independently of 
their properties) closer to the convective envelope; note that Kt 
at a given depth within the star increases only slightly during 
the evolution. Besides, the build up of differential rotation with 
time within the star induces a change in the local Doppler shift 
frequency, which allows a diff'erent damping for waves with dif- 



ferent frequencies and m through the term cr ^ (N^ ~ '^i 



-0.5 



Let us see what these general considerations imply for the 
1.0 Mq model. We start with initial solid body rotation and then 
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Fig. 7. Decomposition of the total flux of angular momentum (solid black) into meridional circulation (long-dashed magenta), shear 
turbulence (dotted blue), and IGW (short-dashed red) in the 1 M© models computed with and without IGW (left and right panels 
respectively). Bold lines indicate negative values for the fluxes Fmc (f), ^s (r), or Fiqw (r), when angular momentum is transported 
towards the central regions by the corresponding mechanism; in the case of meridional circulation and of shear turbulence this 
corresponds respectively to clockwise currents {U2 > 0) and to a positive O gradient. The profiles are shown at three diff'erent ages 
along the PMS. Shaded areas correspond to convective regions 



follow the transport of angular momentum when the radiative 
layers appear At that moment differential rotation has not yet 
developed, and the local frequency cr of individual waves in the 
very thin radiative zone is similar to their emission frequency 
0) at the base of the convective envelope. However slight dif- 
ferential rotation soon builds up as a result of stellar contrac- 
tion along the Hayashi track, which induces a Doppler shift be- 
tween the emission and local IGW spectra. As a consequence, 
low-frequency low-degree waves, which undergo the largest dif- 
ferential damping between retrogade and prograde components, 
soon penetrate all the way to the central regions where they de- 
posit their negative momentum and very efficiently spin down 
the core whose amount of angular momentum is minute (see 
Fig. 6). This explains the strong positive gradient in the pro- 
file of Q. below ~ 0.2 Rq, while the negative gradient of Q in 
the external layers results from ongoing stellar contraction. As a 
consequence a peak builds up in the internal rotation profile with 
a core spinning at lower rate than the stellar surface all along the 
PMS. 

We show in Fig. 7 the total flux of angular momentum car- 
ried by the waves as a function of depth within the 1 Mq model, 
and compare it to the contribution of meridional circulation and 
shear turbulence at different evolution stages. We note first that 
the transport of angular momentum is generally dominated by 
the waves in the radiative layers where they can propagate, ex- 
cept in the early times when meridional circulation dominates 
in the most external regions below the convective envelope (up- 
per panel at 14 Myr). Since downward propagating waves are 
totally damped as soon as the condition Q.{r) = oj/m + ficz is 
fulfilled near the corotation radius, the total flux Fjcw drops and 
remains negligible below the O peak. This can be clearly seen 
in the middle and lower panels in Fig. 7 at 35 and 55 Myrs; at 
that time meridional circulation dominates in the regions below 



~ 0.15 and 0.2 M© respectively, while IGW are dominant in the 
outer regions. Note that the total flux of angular momentum is 
dominated by IGW when they are accounted for and is larger 
by two orders of magnitude compared to the case without IGW. 
Overall, IGW do shape the circulation patterns, leading to the 
appearance of several loops in the whole radiative region as can 
be seen in Figs. 7 and 8. 

Let us add a final remark. As explained in § 2.3.1, we have 
increased by a factor 2 the IGW luminosity in order to account 
for the results by Lecoanet & Quataert (2012) who predict the 
IGW flux due to turbulent convection to be a few to five times 
larger than in previous estimates by e.g. Goldreich & Kumar 
(1990) and Goldreich et al. (1994). In order to test the impact 
of this assumption, we have computed two additional models for 
the 1 M0 rotating star with multiplying factors of 1 and 5. We 
find that this has no impact on the conclusions, as can be seen in 
Fig. 10 where we plot the corresponding rotation profiles at the 
arrival on the ZAMS. 



5.2. Impact of the stellar mass 

For all the stars within the considered mass range, strong dif- 
ferential rotation with a fast rotating core is obtained under the 
combined action of stellar contraction and meridional circulation 
when IGW are not accounted for. Besides, in all cases IGW do 
break-up the stellar core, which results in a peak in Q. at r~0.25- 
0.3R, as in the 1 Mq case. This can be seen in Fig. 10 where 
we show the rotation profiles at the anival on the ZAMS for all 
our models (black and red lines correspond respectively to the 
models computed without or with IGW). 

Let us note however that the impact of IGW is slightly dif- 
ferent in stars more massive than ~ 1.6 M©. This is illustrated 
for the 2 M© star in Fig. 1 1 where we decompose the total flux 
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Fig. 8. Meridional circulation currents in the 1 M© models computed with and without IGW (top and bottom respectively) at three 
different evolution ages along the PMS (14, 35, and 55Myrs from left to right)). Blue and red lines indicate clockwise {U2 > 0) and 
counterclockwise {U2 < 0) circulation respectively. Hatched areas correspond to convective regions 



of angular momentum within the model according to the var- 
ious transport processes at three different ages, and in Fig. 12 
where we follow the corresponding evolution of the radial pro- 
file of Q.. For this more massive star, IGW emitted by the con- 
vective envelope dominate during the first part of the PMS and 
manage to slow down the most central regions as in the 1 M© 
case (top panel. Fig. 11). However those waves fade away when 
the convective envelope becomes too thin and are supplanted 
by those emitted by the convective core at the approach of the 
ZAMS (see Fig. 4). During that transition period (middle panel 
in Fig. 11), meridional circulation dominates the transport of 
angular momentum although shear turbulence also contributes 
more efficiently near the most central regions (between 0.05 and 
0.1 Mq) and in the most external layers; as a result, the core 
slightly accelerates and eventually manages to rotate faster than 
the outer radiative layers, although not fast enough for the peak 
to be erased. Once the convective core has sufficiently developed 
(lower panel. Fig. 11), the IGW emitted in the central regions 
will start conveying angular momentum very efficiently towards 
the core; at that time meridional circulation remains however the 
dominant process in the most external radiative layers. 



6. Global stellar properties, surface rotation and 
lithium abundance 

We summarize in Table 1 the main properties of our models 
computed under various assumptions. We also include the pre- 
dictions for two additional models of 1 Mq that account for the 
hydrostatic effects of rotation (i.e., the effects of centrifugal ac- 
celeration on effective gravity) and show in Fig. 15 all the cor- 
responding evolution tracks for this star. We see that the rotating 
tracks without hydrostatic effects are hardly modified compared 
to the standard case, the main shift to slightly lower effective 
temperature and luminosity (that implies slightly longer PMS 
lifetime) being due to the effects of the centrifugal force and not 
to rotation-induced mixing. This is in agreement with the predic- 
tions by Eggenberger et al. (2012) (see also Pinsonneault et al. 
1989; Martin & Claret 1996; Mendes et al. 1999). However the 
hydrostatic effects are modest and our general conclusions on 
the evolution of the internal rotation profile and on the impact of 
IGW are not affected by this simplification. We can also note in 
Table 1 that the models computed with IGW have longer PMS 
lifetimes than the others. This simply results from the higher to- 
tal diffusion coefficient for chemicals in the deep radiative layers 
close to the convective core when central H-burning sets in close 
to the ZAMS (see Fig. 14). 
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As shown in Fig. 13, the evolution of surface rotation for the 
models with IGW accounts well for the mean rotation rates col- 
lected by Gallet & Bouvier (2013) for PMS stars in young open 
clusters in the considered mass range. The rotation velocity at the 
arrival on the ZAMS is slightly higher (by a few %; see Table 1) 
in this case than in rotating models without IGW, due to the dif- 
ferent efficiency of the redistribution of angular momentum by 
the various transport mechanisms within the star as discussed 
previously. Again, the hydrostatic effects are negligible. 

The surface lithium abundance at the ZAMS is not signif- 
icantly different in the rotating models without and with IGW, 
as can be seen from Table 1. Indeed this quantity mostly de- 
pends, on one hand, on the temperature at the base of the con- 



vective envelope, which is unaffected since the evolution tracks 
almost superpose, and on the other hand, on the diffusion coef- 
ficient Deff (Eq. 15) in the external radiative layers shown in 
Fig. 14. Since the gradient of Q in the outer part of the star 
is dominated by stellar contraction and is very similar in the 
cases with and without IGW (see Figs. 6 and 12), the result- 
ing Li abundance at the ZAMS is unaffected. The rotating mod- 
els including the hydrostatic effects have slightly higher lithium 
abundance on the ZAMS, in agreement with the behavior found 
by Eggenberger et al. (2012). In a future work we will revisit 
PMS Li depletion taking the influence of the disk lifetime, of the 
initial rotation velocity, and of magnetic braking into account . 
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Fig. 13. Evolution of the surface rotation rate as a function of time for the models computed with IGW. Stellar masses are indicated 
on the tracks. In the left panel, the theoretical predictions for the 0.9, 1 , and 1 .2 Mq models are compared with the observed rotational 
distribution for stars with estimated masses between ~ 0.9 and 1.1 Mq in young open clusters from Gallet & Bouvier (2013) 
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Fig. 14. Diffusion coefficients associated to meridional circulation (red), horizontal and vertical turbulence (black and blue re- 
spectively). The total diffusion coefficient for the chemicals (magenta) and thermal diffusivity (cyan) are also shown. The figures 
correspond to the 1 M© models with and without IGW (left and right respectively) at different evolution ages along the PMS. 
Hatched areas indicate the convective regions 



7. Conclusions 

In this paper we have analyzed the transport of angular momen- 
tum during the PMS for solar-metallicity, low-mass stars (with 
masses between 0.6 and 2.0 Mq) through the combined action of 
structural changes, meridional circulation, shear turbulence, and 
internal gravity waves generated by Reynold-stress and buoy- 
ancy in the stellar convective envelope and core (when present). 



For all the stellar masses considered, IGW are efficiently 
generated by the convective envelope with a momentum lumi- 
nosity that peaks around T^ff ~6200 K, as in the case of main 
sequence stars. These waves soon become an efficient agent for 
angular momentum redistribution because they spin down the 
stellar core early on the PMS, while structural changes lead to 
a negative differential rotation in the outer stellar layers as the 
star contracts. On the other hand, IGW generated by the con- 
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Fig. 10. Rotation profile on the ZAMS for all the stellar masses 
between 0.6 and 2.0 Mq in the cases without and with IGW (full 
black and red dotted lines respectively). In the 1 M© panel, the 
blue long-dashed and the green dashed lines correspond to com- 
putations made with multiplication factors for IGW luminosity 
of one and five respectively, all the other models with IGW being 
computed with a multiplication factor of two (see Eq. 9) 
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vective core close to the arrival on the ZAMS carry much less 
energy, except in the case of stars more massive than ~ 1.6 Mq. 
Over the whole considered mass range, IGW were found to sig- 
nificantly modify the internal rotation profile of PMS stars and 
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Fig. 15. Impact of rotation, IGW, and of the hydrostatic effects 
on the evolution track of the 1 Mq star. The square indicates the 
point where the radiative zone appears and internal transport of 
angular momentum starts 



lead to slightly higher surface rotation velocity compared to the 
case where only meridional circulation and shear turbulence are 
accounted for. 

The exploratory results presented in this paper show 
the ability of IGW to efficiently extract angular momen- 
tum in the early phases of stellar evolution, as antic- 
ipated by Talon & Charbonnel (2008) and as shown by 
Charbonnel & Talon (2005) and Talon & Charbonnel (2005) for 
solar-type main sequence stars. We now plan to investigate the 
influence of the disk lifetime, of the initial rotation velocity, and 
of magnetic braking during the PMS over a broader mass domain 
in order to compare model predictions with large data sets that 
are currently being collected to trace the rotational properties of 
young stars. 
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